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Abstract. In this paper we prove the global well-posedness for a three-dimensional Boussi- 
nesq system with axisymmetric initial data. This system couples the Navier-Stokcs equation 
with a transport-diffusion equation governing the temperature. Our result holds uniformly 
with respect to the heat conductivity coefficient k > which may vanish. 

1. Introduction 

The Boussinesq system is widely used to model the dynamics of the ocean or the atmosphere. 
It arises from the density dependent incompressible Navier-Stokes equations by using the so- 
called Boussinesq approximation which consists in neglecting the density dependence in all 
the terms but the one involving the gravity. This system writes 

dtV + V ■ Vv - Av + Vp = pe^, {t, x) G M+ x 
dtp + V ■ Vp = nAp, 
divv = 0, 

V\t=o = Vo, p\t=o = Po- 

Here, the velocity v = {v^,v'^,v^) is a three-component vector field with zero divergence, 
the scalar function p denotes the density or the temperature and p the pressure of the fluid. 
The coefficient k > is a Reynolds number which takes into account the strength of heat 
conductivity. Note that we have assumed that the viscosity coefficient is one, one can always 
reduce the problem to this situation by a change of scale (as soon as the fluid is assumed to 
be viscous) which is not important for global well-posedness issues with data of arbitrary size 
that we shall consider. The term pe^ where = (0, 0, 1)* takes into account the influence 
of the gravity and the stratification on the motion of the fluid. Note that when the initial 
density po is identically zero (or constant) then the above system reduces to the classical 
incompressible Navier-Stokes equation: 

{dtV + V ■ Vf — Av + Vp = 
dhfv = 
V\t=o = Vq. 

From this observation, one cannot expect to have a better theory for the Boussinesq system 
than for the Navier-Stokes equations. The existence of global weak solutions in the energy 
space for ([2]) goes back to J. Leray [28]. However the uniqueness of these solutions is only 
known in space dimension two. It is also well-known that smooth solutions are global in 
dimension two and for higher dimensions when the data are small in some critical spaces; see 
for instance for more detailed discussions. In a similar way, the global well-posedness 
for two-dimensional Boussinesq systems which has recently drawn a lot of attention seems 
to be in a satisfactory state. More precisely global well-posedness has been shown in various 
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function spaces and for different viscosities, we refer for example to [21 [9l [TTl [TTl [HI [TH [201 
[221 [231 [21] . Fo'^ three-dimensional systems few results are known about global existence. We 
can quote the result of R. Danchin and M. Paicu [i8\ who proved a global well-posedness 
result for small initial data belonging to some critical Lorentz spaces. 

Let us recall that it is still not known if smooth solutions with large initial data for the 
Navier-Stokes equations can blow-up in finite time in dimension 3. Only some partial results 
are known. For example, in a recent series of papers [I31[I11IISJ global existence in dimension 
three is established for initial data which are not small in any critical space but which have 
some special structure (oscillations or slow variations in one direction). Another interesting 
case of global existence for corresponding to large initial data but with special structure 
is the more classical case of axisymmetric solutions without swirl. Our aim in this paper is to 
establish the corresponding global well-posedness result for the three-dimensional Boussinesq 
system. 

Before stating our main result, let us describe the classical result for the Navier-Stokes 
equation. It is well-known that the control of the vorticity uj which is the vector defined by 
uj = curlf and solving the vorticity equation 

dtuj + V ■ Vu; — Au; = u ■ Vf 

is crucial in order to get global well-posedness results. According to the classical Beale-Kato- 
Majda criterion |6j the control of the vorticity in L°°) is sufficient to get the global 

existence of smooth solutions. The main difficulty arising in dimension three is the lack of 
information about the influence of the vortex-stretching term u ■ Vv on the motion of the 
fluid. Let us now consider a vector field v which is axisymmetric without swirl, this means 
that it has the form: 

(3) v{t,x) = v''{t,r, z)er + v''{t,r, z)ez, x={xi,X2,z), r = {xl + xl)^, 

where (e^, eg, Cz) is the local orthonormal basis of corresponding to cylindrical coordinates. 
Note that we assume that the velocity is invariant by rotation around the vertical axis 
(axisymmetric flow) and that the component of v about eg identically vanishes (without 
swirl). For these flows, we have: 

UJ = {dzV^ — drV^)eg := OUgCg, UJ ■ Vf = UJ. 

r 

In particular ug satisfles the equation 

(4) dtUJg + V ■ Vujg - Aug + ^ = -^^e- 

The crucial fact is then that the quantity C ■= v" evolves according to the equation 

dt( + vV(-{A + \)( = 
from which we get that for all p G [1 , C)o] 

||C(t)l|L.<||Co||L.. 

It was shown by M. Ukhoviskii and V. Yudovich [30] and independently by O. A. Ladyzhen- 
skaya [25] that these new a priori estimates are strong enough to prevent the formation of 
singularities in flnite time for axisymmetric flows without swirl: the system ([2D has a unique 
global solution for vq G such that ujq,^ G fl L°°. We point out that the result is 
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uniform with respect to vanishing viscosity and thus there is no blowup even for the Euler 
equation. Note that in term of Sobolev regularities these assumptions are satisfied when 
vq E with s > |. This regularity assumption is not optimal and has been weakened in 

[27] for vq G and more recently in [1] for vq G H^. 

Our aim here is to extend these classical results for the Navier-Stokes equation to Boussinesq 
systems. The equation for ^ = cj^/r becomes 

(5) + ■ VC - (Ar + \)C = 

and thus the difficulty is to use some a priori estimates on p to control the term in the 
right-hand side of The rough idea is that on the axis r = the singularity ^ scales as 
a derivative and hence that the forcing term drp/r can be thought as a Laplacian of p and 
thus one may try to use smoothing effects to control it. Nevetheless, when k, = 0, since there 
is no smoothing effects on p, one can hope to compensate the loss of derivatives in the right 
hand-side of (jS]) only by using the full smoothing effects of the heat type equation in the left 
hand side. Note that this kind of estimates does not follow from energy estimates and hence 
the problem that one has to face is that the convection term that has to be handled in the 
process is not negligible: this approach naturally leads to some restriction on the size of the 
data. 

In a global existence result for the system ([1]), with k = 0, was established but under 
some restrictive conditions on the support of the initial density namely that it does not 
intersect the axis r = 0. More precisely, 

Theorem 1.1. LetvQ G be an axisymmetric divergence free vector field without swirl and 
such that ^ E L?. Let Pq E L'^ H L°° axisymmetric and such that supp po does not intersect 
the axis (Oz) and Ilz{supp po) is a compact set. Then the system ([1]), with k = 0, has a 
unique global solution {v, p) such that 

-EL^M+;L'), pEL^M+;L'r}L^). 

r 

Here Hz denotes the orthogonal projector over (Oz). 

Since to bound the quantity ||i^/''"||l^l2 ,-12.2^1 one needs to estimate ||p/r(t)||2,2. The idea 
of the proof was to get an estimate from below for the distance of the support of p to the 
vertical axis (note that this distance remains positive as long as the solution remains smooth 
because of the assumption on the initial density). Note that for this approach, it is crucial 
to have a transport equation for the density, it fails for k > because p cannot be supported 
away from the axis even if the initial data is. 

In this paper, by using a different approach, which uses more deeply the structure of the 
coupling between the two equations of ([I]), we remove the assumption on the support of the 
density and we give a global well-posedness result with uniform bounds with respect to the 
heat conductivity k. Our main result reads as follows. 

Theorem 1.2. Consider the Boussinesq system ([1]) for n > 0. Let vq E be an axisym- 
metric divergence free vector field without swirl such that ^ E L"^ and let po E L'^ n B^-^ an 
axisymmetric function. Then there is a unique global solution {v, p) such that 

V E C(M+; H') n LL(M+; H^) H 5^ J, 
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Moreover the estimates in the above spaces are uniform for k> in any bounded set o/M+. 

The definition of the Besov spaces ^ is recalled below. 
Let us give a few comments about our result. 

Remark 1.3. We find at most an exponential growth for the velocity: there exists Cq > 
depending only on the data such that for every k G [0, 1], we have 

Remark 1.4. The uniformity of the norms with respect to the conductivity k can be estab- 
lished for every k G [0, A\ for every A > 0. Nevertheless, as we shall see below, we need to 
use a different transformation of the equation when n is close to one. 

Remark 1.5. The result of the theorem remains true if we change the Besov space B^^ for 
the Lebesgue space for m > 3. We have not tried to get the best result in terms of the 
regularity of the velocity. At the price of more technicality, it is probably possible as in \\\ 
to get the same result by assuming only that Vq G -ffa. 

Remark 1.6. By using the control of Vf in B^ which is given by Theorem II. 21 on can easily 
propagate by classical arguments higher order regularity for example higher Sobolev 
regularity. 

Let us explain our strategy for the proof. The crucial part in the proof consists now in finding 
suitable a priori estimates for ((^, p). Let us describe the idea in the case that k = (this is 
the one that one needs to understand in priority because of the lack of smoothing effect). As 
we have already noticed the coupling between the two equations does not make the original 
Boussinesq system ([1]) well-suited for a priori estimates. Since the right hand side of ([5]) 
behaves roughly as a Laplacian, we need to fully use the left hand side to control it. The 
main idea is to use an approach related to the one used for the study of two-dimensional 
systems with a critical dissipation, see [221 [21] • It consists in diagonalizing the linear part of 
the system satisfied by and p. We introduce a new unknown F which formally reads 

r = C- A + -a. 'J^p:=^-Cp 

and we study the system satisfied by (F, p) which is given by: 

dtr + vVr-{A+ '^dr)T = -[C,vVp], dtp + vVp = 

where [C, f ■ V] is the commutator defined by 

[C,vV]f = C{vVf)-vVCf. 

We can thus get a priori estimates for F and p (note that they are obvious if we negect the 
commutator) and then use them to deduce estimates for (. The main difficulties that one 
has to deal with are twofold. The first one is the study of the operator C, to make this 
argument rigorous, we need to prove that £ is well-defined and is a bounded operator on 
L^. The second one is the study of the commutator. Once estimates on ( are obtained, the 
situation gets close to the standard axisymmetric Navier-Stokes equations. When k is non 
zero, it turns out that the same kind of transformation can be used and that it depends 
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smoothly on k. The new unknown T ^ = {1 — k)i^ — Cp solves the same convection diffusion 
equation as F. This is due to a surprising commutation property between A and C which 
is stated in Lemma 13.41 Consequently, we are again able to obtain and estimate for C, from 
an estimate for as soon as k 7^ 1. The case that k is close to one (which is easier since 
there is a non vanishing smoothing effect in the density equation) can be handled by using 
a different transformation. 

The paper is organized as follows. In section [2] we fix the notations, give the definitions of the 
functional spaces that we shall use and state some useful inequalities. Next, in section [3l we 
study the operator £, this amounts to study an elliptic equation with singular coefficients. 
In section m we obtain a priori estimates for sufficiently smooth solutions of ([1]), they are 
obtained by using the procedure that we have just described. Finally, in section [5], we give 
the proof of Theorem 1 1.2t we obtain the existence part by using the a apriori estimates and 
an approximation argument and then we prove the uniqueness part. 



Throughout this paper, C stands for some real positive constant which may be different in 
each occurrence and Cq for a positive constant depending on the initial data. Moreover, 
both are assumed to be independent of k for k > in a bounded set. We shall sometimes 
alternatively use the notation X <Y for an inequality of the type X < CY . 
For s G M, we denote by if'^(M^) the standard Sobolev spaces: u belongs to if u is a 
tempered distribution and 



Now to introduce Besov spaces which are a generalization of Sobolev spaces we need to recall 
the dyadic decomposition of the whole space (see [12j). 

Proposition 2.1. There exist two positive radial functions x ^ "DiM?) and ip G T'(M'^\{0}) 
such that 



(2) supp ip{2-P-) n supp ¥?(2-«-) = 0, if |p - g| > 2, 

(3) g > 1 ^ suppx n supp v^(2^'^) = 0. 

For every u G iS'(M^) we define the nonhomogeneous Littlewood-Paley operators by. 



2. Preliminaries 




We shall also use the homogeneous version if^(M'^): for s < 3/2, m G if* if u G L]^^ and 




(1) x(0 + E ^(2-^0 = 1, 3 < X'(0 + E V''(2-«0 < 1 G M3, 




A_iM = x(D)m; Vg G N, AgM = y?(2-''D)M and 





-i<i<g-i 



One can easily prove that for every tempered distribution u, we have 
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In the sequel we will frequently use Bernstein inequalities (see for example |12j). 

Lemma 2.2. There exists a constant C such that for kEN,l<a<b and u G L"- , we have 



and for g G N 



sup II^^VIU" < C^2'i^^+^^-a'-b^^\\Squ\\L^, 

\a\=k 



C-^2''%\u\\La < sup ||9"A,M||ia < C^'2«*-'||AgM||ia. 

\a\=k 



Let (p, r) G [1, +00]^ and s G M, then the Besov space 5^ is the set of tempered distributions 
u such that 



:= (2''1|A,w|U.)^^ < +00. 



We remark that the Sobolev space if* agrees with the Besov space -B|2- Also, a straightfor- 
ward consequence of the Bernstein inequalities is the following continous embedding: 

(7) BI^^^^bX[^-~^'\ Pi<P2 and n < ra. 

For any Banach space X with norm || • \\x and fonctions f{t,x) such that for every t, 
fit,-) G X, we shall use the notation = II ||/||x||lp{[o,t])- 

A useful application of Besov spaces is the following logarithmic estimate for convection 
diffusion equations. 

Proposition 2.3. There exists C > such that for every n > 0, p E [1, 00] and for every p 
solution of 

(dt + v-V - kA)p = f, p(0, x) = po{x) 
with V a divergence free vector field, the following estimate holds true 

||p(t)||BO^<C(||po||BO, + ||/||LiBO,)(l+ / \\Vv{T)U^dTy Vt>0. 



We refer to [21] for the proof. Note that the amplification factor is only linear in ||Vf ||loo. 

3. About an elliptic problem 

The aim of this section is the study of the operator £ = (A + |)^^^. This is the heart of 
the paper since this is crucial to make rigorous the argument sketched in the introduction. 
This amounts to study the regularity of the solution of an elliptic equation with singular 
coefficients. This is the goal of the following proposition. 

Proposition 3.1. Let p G H'^{M.^) axisymmetric, then there exists a unique axisymmetric 
solution f G of the elliptic problem 



Moreover, for every p G [2, +00), there exists an absolute constant Cp > such that: 

(9) II/IIl. <Cp||p|Up. 



The important fact in this proposition is the U' estimate ([9]) which only involves the norm 
of p. An immediate consequence of Proposition 13.11 is that C defines a bounded operator on 
for every p G [2, +cx)). The additional regularity is used to give a meaning to the 
equation. Because of the 1/r singularity on the axis, we cannot give a meaning to the term 
drf /r as a distribution when / is merely V . When / is in H^, there is no problem we have 
that drf /r G L]^^ since for every compact set C M'^ 

ll*^^.^!! ^ IIV7/II ll^ll 

\\ — UiK)< llV/lUe 11-11^6^^^ <+00 

thanks to the Sobolev embedding C in dimension 3. 



Proof. Let us first prove the existence of a solution satisfying the required properties. We 
can first assume that p G C^(M^) and then conclude by density. Since the elliptic operator 
has some singular coefficients, we shall use an approximation argument. Since we have by 
definition that rdr = x/^ ■ V with the notation Xh = (a;i,X2,0), we shall consider for e > 
the elliptic problem 

(10) (^ + ^^^-v)/ = ;:^x.-Vp. 

Since the coefficients are not singular any more, there is a unique solution for this problem 
given by the classical methods. By standard regularity arguments, this solution is in the 
Schwartz class and hence the following a priori estimates are justified. Moreover, since p is 
axisymmetric, is also axisymmetric. 

We shall first prove that the solution of flTUl) satisfies the estimate ([HD with a constant 
independent of e. In the proof of the a priori estimate, we shall denote by / for notational 
convenience. By taking the L^(]R^) scalar product of (ITO!) with (r^ + e)|/|^^^sign(/), we find 

(11) / Af\frhign{f){r' + e)dx + 2 [ (x, ■ V/) |/r^sign(/) da; 
= [ Xh-Vp\frhign{f)dx. 

For the first term in the left-hand side, an integration by parts yields 

/ Af\fr'signif){r' + e)dx = -{p - I) [ \fr'\Vff{r' + e)dx 

-2 / \frhign{f)xh-Vfdx. 

Consequently, we get that 

(12) / Af\fr'sign{f)dx + 2 [ {xh-Vf)\fr'sign{f)dx 
= -ip-l) [ \fr'\\/fmr' + e)dx. 



For the right-hand side of (fTTI) . we also obtain from an integration by parts that 
/ Xf,-Vp\frhign{f)dx = -2 f p\frhign{f)dx 



-ip-1) [ pxH-vf \fr'dx. 



By using the Holder inequality, this yields 



Xh-yp\f\'' ^sign(/)rfx 



< 2 IIpI 



LP 



p-1 
LP 



2-2 

2 II /ill 

LP WPWlp 



\Vnf?\fr^T^dx 



where Vh = (c^i,f^2)*- By using this last estimate, ffTTl) and ffT^ . we obtain that for every 
e>0 



/ \fr'\vf\'r'dx < 2||piu.||/r^;^ 



p-2 



Consequently, by using the Young inequality, we get that 
(13) [ \fr'\Vf\'r'dx<c(\\ph. 



\^hf\'\fr'r'dx 



LP ^ WPWlp 



p-2 
LP 



for some C > independent of e. To conclude, we can use the following inequality: for 
p G [2, oo[, we have 

(14) 



This is a special case of Caffarelli-Kohn-Nirenberg inequality [lOj (we shall recall the proof 
will be given in the end of the proof of the proposition). From f[T51) and ([SD, we obain that 



<c{\\p\\lp II/I|l. + ||pIIL) 



where C is independent of e and and thus we obtain the estimate for the solution of ffTUl) 
by using the Young inequality. 

When e goes to zero, we get the existence of / G and a subsequence e„ such that the 
solution of (fTOj) converges weakly to / which satisfies the desired estimate 



LP 



< Cp\\p\ 



LP- 



In order to get that / solves the equation ([H]), we need more information on Indeed, the 
difficulty is to give a meaning to the term d^f /r which is not well-defined as a distribution 
when / is merely U'. We shall use a uniform estimate for the solution of ffTOl) but which 
also involves the norm of p. This is why we have required more regularity on p. At first, 
let us notice that since we assume that p is in L^, we have that is uniformly bounded in 
L^. Next, we multiply ffTOl) by A/, we get that 



(15) 



\/\f\^dx + 2 I ^I^A/dx 



r"^ + e 



Afdx. 



The right-hand side can be estimated by 

Xh ■ Vp 



+ e 



Afdx 



< \\^''P\\ IIAfll 

< II — \\m\'^i\\L-^- 



r 



Since we assume that p is axisymmetric, we can use that 

12^2 , ^ 



Ap - - dip = ^dlp + ^dlp - 2^a?,p 



r 



and hence that we have the estimate 



drp I 



^2 < 4||p||h2- 



This yields 
(16) 



Xh ■ Vp 



< 4||p||^2||A/||i2. 



Next, we can study the second term in the left-hand side of (|T5|) . We have by integration 
by parts that 



+ e 



Afdx 



Xh 



+ e 



V(|V/n-2/ (V(-^).V)/.V/ 



/ V-(^)|V/r-2/ (v(^)-v)/-V/rfx. 



Next, we infer 



V2 + e^ 



2e 



and since / is axisymmetric 

(v(;5^)-v)/.V/^(-J 

This yields 

Xh-Vf 



2r' 



+ e (r2 + £)2 



\drf\' 



(17) 



■A/ rfx = 2 



:\drf\'dx>0. 



Consequently, we get from (iTSl) . f[T6|) and ffTTI) that 

IIA/IU2 <4||p||^,2 

and hence that 

(18) ||/||i/2<C||p||H2 

with C independent of e. From this uniform H"^ estimate for we get that (up to a 
subsequence not relabelled) converges weakly in to some / G and then that / is a 
weak solution of ([H]) • This ends the proof of the existence of a solution. The uniqueness is a 
consequence of the standard energy inequality. 
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Let us now come back to the proof of (fT^ . Since div{xh) = 2, by integrating by parts, we 
obtain 



= -P [ irdrf)\frhignif)dx 
Jr» 

< p( [ IrdrfWr'dx 



E 
2 

LP- 



Therefore we find the desired estimate, 



Ip<j I IrdjWr'dx. 



p 

This ends the proof of Proposition 13.11 

□ 

In the proof of the main result, we shall also need to use the operator (A + ^dr) ^dz/r. The 
aim of the following proposition is to define rigorously this operator. 

Proposition 3.2. Let p G L^(M^) be axisymmetric such that dzp/r G L^(R^), then there 
exists a unique axisymmetric solution f G of the elliptic problem 

(19) (A + h)f = ^. 

Moreover, for every p G [2, +oo), there exists an absolute constant Cp > such that: 

(20) II/IIl. <Cp||p|U.. 

Again, the important fact is the estimate (I2UI) which only involves the norm of p. From 
this estimate, we get that the operator (A + ^dr) dz/r is a bounded operator on L^. The 
additional regularity p G L^, dzp/r G is again only used to get the regularity on the 
solution which allows to give a meaning to the equation. 

Note that the assumptions on p here are different from the one of Proposition 13. 1[ This 
comes from the fact that we shall need to use (A + ^(9^.) when p is a smooth solution of 

([T]) whereas, we shall only need to use (A + ^5^.) (f ''p) where (f , p) is a smooth solution 
of (lH). In the first case, as we have seen in the proof, the regularity on p ensures that drp/r 
is in L^. In a similar way, in the second case for a smooth solution of ([T]) such that uj/r G L^, 
we indeed have that ^{v^p) G L^. 

Proof. The proof follows the same lines as the proof of Proposition 13. 1[ Consequently, we 
shall just indicate the main difference. We now consider the regularized problem 

(21) (A + ^x.-V.)/ 



r'^ + e yV2 + e 

By multiplying the equation with (r^ + e)|/p~^sign(/), we get by using again (fT2l) and an 
integration by parts for the right-hand side that 



/ \fr'\^f\\r' + e)dx< [ V^:^\p\\fr'\dzf\ 



dx. 
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From the Holder inequality, we obtain that 

( / i/r'iv/ivdx)^<iipiu.ii/iig' 

and hence by using again (fn|) . we finally get that 

II/IIl. < CpIIpIUp. 

This will give the estimate ( l20i) by passing to the limit. 

To prove an estimate on /, we again multiply ( 1211) by A/, since / is axisymmetric, we 
get from flTTI) that 

II A f II ^ ll^^^ll 
r 

and this provides the estimate for /. We can then pass to the limit to get a solution of 
f ll9p with the claimed properties. The uniqueness follows from the classical energy estimate. 
This ends the proof of Proposition 13.21 

□ 

The aim of the next two lemmas is to prove some identities involving the operator C = 
(A + ^dr) ^ which will be useful to get the equation satisfied by Cp in order to diagonalize 
the system. 
At first, we have 

Lemma 3.3. For any smooth axisymmetric function f we have the identity : 

Proof. We first obtain that 

1 11 

-drrf = dr{-drf) + ^9^/ 

= dr{dr{^ + ^)+\drf 

^ ly ry ^ ' ry ^ 

= 5rr(-) + 4<9^/ - 4/ 

ly ly ^ ryO 

= (a„. + 'a.)({) 

\ ty ' ry ry ry ry 

It follows that 

\ ry ' ry ry ^ ry ' ry ry ' ry ' ry 

ry ry ^ ry ' ry 

Note that we have used the fact that the operators (A + jdr^ ^ and dz commute since the 
coefficients of the operator A + "^d^ do not depend on the variable z. This is the desired 
identity and therefore, this ends the proof of Lemma [3. 3[ □ 

We shall also use the following: 
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Lemma 3.4. For every smooth axisymmetric function p, we have the identity: 

£ Ap = ( A + -dr) Cp. 
Proof. At first, by direct computations, we find that 

(22) [A, -dr] = -2{\d', - ^dr) = --dr{-dr • )• 

Now, let us consider f = C p. By definition, / solves the elliptic equation 

{A + -dr)f=-drp. 

Consequently, we get that m = (A + ^dr)f solves the elliptic equation: 
By using the formula (1221) . we obtain for the right-hand side 

(A + -dr) (-drp) = -dr(Ap + -drp) + [A, -dr]p = -9, (Ap + -9,p - -drp) = -drAp. 

\ ly* ' ^ /y ' /y ^ /y ' <- ry ry ^ ry ry ' ry 

This proves that u solves the equation 

2 1 
(A + -dr)u = -drAp 

^ ry ' ry 

and hence that u = CAp. Since m = (A + ^dr)f = (A + ^dr)Cp, this ends the proof of 
Lemma 13.41 

□ 

4. A PRIORI ESTIMATES 

This section is devoted to the a priori estimates needed for the proof of Theorem 11.21 We 
shall prove two results: the first one deals with some basic energy estimates. The second 
one which is more difficult deals with the control of some stronger norms. 

Proposition 4.1. Let (f,p) be a smooth solution of ([1]) then 

(1) for p E [1, oo] and t G M+, we have 

\\p{t)\\Lp < IIpoIIlp, 

(2) for Vq G L^, Po G and t G M+ we have 

\Hml^+ f\\Vvir)\\Ur<Co{l + t'), 
Jo 

where Cq depends only on ||fo||L2 and ||po||l2- 

Note that the axisymmetric assumption is not needed in this proposition 

Proof. The first estimate is classical for convection diffusion equations with a divergence free 
vector field. For p = oo it is just the maximum principle, while for finite p, it is a consequence 
of the fact that 

{dt + v - V - kA)\p\^ <0. 
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For the second one we take the L^(M^) scalar product of the velocity equation with v. From 
integration by parts and the fact that v is divergence free, we get 

(23) ^lii^^'^^ii^^ + < Mt)\\L4pmL^- 

This yields 

jtMt)h^ < IIpWIU- 

By integration in time, we find that 

\\v{t)\\L2 < \\vo\\l2 + [ \\p{T)\\L2dT. 

Jo 

Since ||p(t)||L2 < ||po||l2, we infer 

||^^(t)||L2 < ll^^olU^ +i||Po|U2- 

Plugging this estimate into ( j23l) gives 

1 
2 

This gives the second claimed estimate and ends the proof of the proposition. □ 
We shall next prove the following result. 

Proposition 4.2. Let vq G H^, with ooo/r G and po ^ flL^. Then any smooth solution 
(f , p) of ([T]) with p axisymmetric and v axisymmetric without swirl satisfies: 

(1) for every t G R+ 



- \\v 



+ j ||Vt;(r)||^2rfr < -||wo||i2 + (||t;o||L2 +t||po||L2)||po||L2^• 



-it) 

r 



L2 



(2) for every t G M+ 

11^(^)11^1+ f \Hr)\\l2dr < Coe''^^ 
Jo 

where Cq depends only on the norms of the initial data. 

Note that the axi-symmetry is crucial in this proposition. The estimates are uniform for 
K > in a bounded set. 

Proof. We shall use the notation ( = ug/r where the vorticity u is given by u; = ugeg since 
the flow is axisymmetric. The equation for ( reads 

(24) ^d, + vV)C-{A + \)C = -^- 
By using the operator £, this equation can be written as 

(25) {dt + v-V)C-{A + ^9,) (C - Cp) = 0. 

Applying C to the equation for p we first get 

dtCp + V ■ VCp — kCAp = —[C,v- V]p 
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and hence by using Lemma 13.41 we find 

2 

(26) dtCp + v- VCp - k(A + -dr)Cp = -[C,v ■ V]p. 



In view of fl25l) and fl26l) . we can set 

r := (1 - k)C - Cp. 

We find that F solves the equation 

{dt + v- V)r - (A + '^dr)T = div (vCp) - C{v ■ Vp). 

Taking the L^(M^) inner product with F and integrating by parts in the usual way we get 
since v is divergence free that 

(27) ^|l|r(t)|U. + l|vr(t)||i. < ||vr(t)|U.|l^^p||L^ 

- / C{v-Vp)Tdx 
= I + 11. 

To estimate the first term we use successively the Holder inequality, Proposition 13.11 and the 
first estimate of Proposition 14.11 to get 

Ib^CpllLS < ||^||l6||>Cp||l3 < ||^;||i6||p||J:3 < ||t;||i6||po|U3 

Now from the Young inequality, we get 

(28) I<^||Vr||i.+C||t;||i,||po||i3. 

Towards the estimate of the second term in the right-hand side of fl27j) . we can use since v 
is divergence free that 

vVp = dr{v''p) H p + i9^(w». 

r 

Thus we obtain 

V 

C{vVp) = Cdr{v''p) + cC—p) 
+ d,C{v'p). 



Next, we can use Lemma [3.31 to get 

+ d,C{v'p). 



C{vVp) = -p-9,(A + -a,)"'(\(t;») 



Thus we find that 



II = - [ —pTdx+ I (A + -9,) \-d,{v^p))dj:dx 

j ^ j ry ' ry 



C{v'p)d,Tdx 
J1 + J2 + h- 
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From Cauchy-Schwarz inequality, we first obtain that 

IJ2 + J3I < ll^.riU^f (A + -dry\-){v'-p) +\\Civ'p)\\L2) < \\VT\\l2{\\v'-p\\l2 + \\v'p\\l2) 
\ r r L'^ / 

where the last estimate comes from Proposition 13.21 and Proposition 13.11 
Next, by using successively the Holder inequality, the Sobolev inequality 

(29) ll/lUe < IIV/IU. 

in dimension 3, and the first estimate of Proposition 14. 1^ we infer 

\J2 + J3\ < l|Vr||i2||i;||j^6||p||i3 

< ||vr||i2||Vf ||l2||poIIl3. 

In a similar way, we estimate Ji: 

\Ji\ < llrpll 6||i;7r|U6 

< ||r||z,2||p||i3||f''/r||i6 

< l|r||L2||po|U3||t;7^IU6. 

To estimate ||f'"/r||i6, we can use the inequality 

\v"/r\ < j^-^{\uJe/r\) 

which is very useful in the field of axisymmetric solution of incompressible fluid mechanics 
equations, we refer to f2^ ITB] for the proof. Next, since G L5'°°(]R3~)^ gg^ from the 
classical Hardy-Littlewood-Sobolev inequality that 

(30) Wv'/rh^ < Wooe/rh^ = IICIU^. 
From the definition of F we have that 

(31) \\C\\L^<\K-l\-'[\\Th2 + \\CpU.). 

Note that this estimate is uniform for k > and far from k = 1. We will see later how to 
obtain uniform estimates around the value k = 1. Therefore, by using again Proposition 13.11 
and the first estimate of Proposition 14.11 we obtain 

ll^V^lUe ^ l|r|U2 + \\^p\\l2 < ||r||L2 + \\p\\l2 < \\t\\l2 + Wpoh^- 

Consequently 

l-^il ^ l|r||L2||pollL3(||F||i2 + ||poIIl2)- 

Combining theses estimates with the Young inequality we get 

^ l|VF||i2||Vt;||L2||po||i3 + ||F||i2||po||L3(||r||L2 + ||po||l2) 

(32) < i||VF||i2+C||V^;||i2||po||i3 + C||po||L3||F||i2 + C||F|U2||po|U3||po|U2. 



It follows from (|27|), ([28]) and (|32|) that 

^I|r(t)||i2 + ||VF(t)||i2 < ||Vt;||i2||po||i3 + ||po||L3||F||i2 + ||F|U2||po|U3||po||l2 

^ IIpo||l3||F||^2 + ||Vt;||i2||po|li3 + ||po||L3||po|li2. 
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We can then integrate in time and use the energy inequahty of Proposition 14. II - (2) to get 
I|r(t)||i2+/ ||Vr(r)||i.rfr < ||po|lia||Vt;||i.^. + ||po||L3||po|li2t 

+ ||po||l3 [ \\T{T)\\l.dT 



Jo 

By using the Gronwall inequahty we find 

l|r(t)||i.+ r||Vr(r)||i.c/r<Coe^«^ 

^0 

It follows from a new use of flMl) that 

(33) ||CWIlL^<Coe'^°*. 
This proves (1) in Proposition 14.21 

To prove (2), we can now perform an energy estimate on the equation (j4]) satisfied by ug. 
By taking the L^(]R^) scalar product of (jl]) with ug we get 

T^^W^eit^h + WVuJeWh + \\iue/r\\'l2 < f v''{uje/r)ujedx- [ drpujedx. 

Thanks to an integration by parts, we have that 



drpoJedx 



< ||p||l2 (||Vu;e||L2 + \\ujg/r\\L2) 



and hence, by using the Holder inequality and the Sobolev inequality (l29l) . we find 

^^ll^eWlli^ + \\^uje\\l2 + \\uJe/r\\l2 

< \\v\\L3\\uJe/r\\L2\\^e\\L(i + IIpIIl^ (|| Vt^elUa + \\ujg/r\\L2) 

< \\v\\L3\\ujg/r\\L2\\Vuj0\\L2 + ||p|| L2 (|| Vc^e || L2 + H^^e/r || l2 ) . 

Thus from (133|) . the first estimate of Proposition 14. II and the Young inequality, we obtain 

(34) ^II^^WII'^ + I{\\'^^e\\l2 + \\ujg/r\\l2) < UpoHi^ + Coe^'"\\vit)\\l,. 

By interpolation, the Sobolev embedding (129|) and the second estimate of Proposition 14.11 
we have 



f\\v\\L4'^vh.<Co{l+f). 

Jo 



La < / ||i;||l2||^^||l6 < 
Jo 

Therefore we get by integrating fl5^ in time that 



K(^)lli2+ / (||Vc.,(r)||i. + ||5^(r)||i.)^r<e^ 
Jo ^ 



Since we have ||l(j||l2 = ||ci;6i||l2 and 

I|Vu;|Il2 = ||Vu;e||i2 + \\uje/r\\\2, 
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we finally obtain that 

\Ht)\\h+ f\\Vu{r)\\UT<e^^\ 
Jo 

This ends the proof of Proposition l4.2l where k belongs to a compact set that does not contain 
1. Let us now see how to get uniform bounds around k, = 1 which is more easy and does not 
require the use of the operator C We write the equation of the density under the form 

2 2 
dtp + V ■ Vp — (A + -dr)p = (1 — K)Ap drp. 



We set Ti = ( — ^. Then combining this equation with fl24j) we find 

Taking the scalar product of this equation with Fi and integrating by parts we get 

< ^(^nivp||i. + ^||vr,||i.. 

Integrating in time yields 

l|ri(t)||i. + ||vri||i.,, <c(«:-i)2||Vp||i.,.. 

Combining this estimate with the energy estimate 

'^^I|Vp||l2^2 < ||po||l2 

gives 

l|ri(t)||i. + ||vri||i.,,<c^^^||polli.. 

This gives the desired result. □ 



The following proposition gives some more precise information than stated in Theorem 11.21 
about the solution. This will be useful to prove the uniqueness result. 

Proposition 4.3. Let vq G be a divergence free axisymmetric without swirl vector field 
such that uoo/r G and po & (1 L"^,m > 3 or po G fl B^^ an axisymmetric function. 
Then any smooth solution {v, p) of the system ([1]) satisfies for every p g]3, oo] 

1+3 + ||Vt;||riL- < Coe^"'. 
The estimate is uniform with respect to k, lying in a bounded set. 

Proof. We first prove the result in the case of po E L"^ H ^. Let g G N and set Vg := AgV. 
Then applying the operator Ag from the Littlewood-Paley decomposition to the velocity 
equation and using Duhamel formula we get 

(35) vg[t) = e'\{0) + [ e^'-^'^^ AgV{v ■Vv){T,x)dT + [ e(*-")^A,P(pe,)(r, x)cir, 

Jo Jo 

where V is the Leray projection on divergence free vector fields. Now we will use two 

estimates: the first one is proved in [12] 
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The second estimate is 

||A,P/|U. <C||AJ|U.. 

This last estimate is a consequence of the fact that AqV = ip{2~'^D) with ip G P(]R^). 
Therefore, we get from (!35|) that 

\\vg\\LiLp<2-''\\v,{0)\\L^ + 2'"' f \\A,{vVv){T)\\L.dr + 2-"^\\A,p\\LiL,. 

Jo 

It follows from the above inequality, the Besov embeddings ([7]) and Proposition 14.11 that 

ll^^ll ^ 1+1 < \\A^iv\\lilp + \\vo\\ -1+1 + 11^^- Vt;|| + IIpII _i+| 

< Ct\\v\\L^L^ + \\vo\\hi + \\V'S)V\\ 3 + IIpII rifiO 

< Co{l + t^) + \\v®v\\^^3 +\\p\\l}bO- 



3 

Since B21 is an algebra, we have 



If ® f II 3 < C||f II 3 • 



Moreover, the embedding ^ B2 i combined with the second estimate of Proposition 14.21 
gives 



'L^Bl, 



Consequently, we obtain 

(36) \\v\\ <Coe^''*+ii/^iiLiB30,- 

r 1 p P t d.l 

It remains to estimate the norm of the density. For this purpose, we first use the logarithmic 
estimate described in Proposition 12.31 and ([7j) to get 

\\p{t)Ul^ < C||po||BO,(l + £ll^(r)||B^/r) 
< q|po||BO,(l+ flHrn i+|rfr). 

, 1 

Set V(t) := llf II 3 , then combining this estimate with (l36l) yields 

V{t) < Coe^-^' + llpollijo fv{T)dT. 
We conclude now by Gronwall lemma. 

Let us now come back to the case that po E L"^ (1 L™, with m > 3 which is more easy than 
the previous case. The same proof as above for the velocity yields 

ll^^ll i+a < C'oe^"* + IIpII _i+3 
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and it still remains to estimate the density. Let us set mi := min(m,p) > 3 then by Besov 
embeddings and the first estimate of Proposition 14.11 we get 

IIpII -1+1 ~ IIpII -1+7^ ~ IIpIIlil™! ^^IIpoI|l'"i < i||pollL2nL™- 

Note that the last estimate holds by interpolation. This ends the proof. □ 

5. Proof of the main result 
For the existence part of Theorem 11.21 we smooth out the initial data as follows 

where Sn is the cut-off in frequency defined in section 2. We start with the following stability 
results. 

Lemma 5.1. Let vq be a free divergence axisymmetric vector-field without swirl and po (^n 
axisymmetric scalar function. Then 

(1) for every n G N, fo,n CLnd po^n are axisymmetric and div i;o,n = 0. 

(2) If Vq G is such that (curl VQ)/r G L"^ and po ^ L"^ Cl B^-^^. Then there exists a 
constant C independent of n such that 

WvqAh^ < ll^oll/fi, ||(curl t;o,„)/r||^2 < C||(curl VQ)/r\^^^, 

||P0,n||L2 < ||Po||l2, llPo.nllsO^ < C'llPoHfiO^- 

Proof. We have vo,n = 2^"'x(2"-) * "^0. The fact that the vector field fo,n is axisymmetric is 
due to the radial property of the functions for more details see [4]. The estimate of fo,n 
in is easy to obtain by using the classical properties of the convolution laws. The proof 
of the second estimate for the velocity is more subtle and we refer to [T] where it is proven 
in the general framework of Lebesgue spaces, that is in L^, for all p G [1, cxd]. The estimates 
for the density follow by standard convolution inequalities. □ 

We have just seen in Lemma 15.11 that the initial structure of axisymmetry is preserved for 
every n and the involved norms are uniformly controlled with respect to this parameter n. 
Thus we can construct locally in time a unique solution (f„,Pn)- This solution is globally 
defined since the Lipschitz norm of the velocity does not blow up in finite time as it was 
stated in Proposition 14.31 Note that since po ^ B^i (or ""^ > 3), the norm of po,n 
is also uniformly bounded thanks to the continuous embedding B*^-^ C L^. By standard 
arguments we can show that this family {vn, pn) converges to (f,p) which satisfies in turn 
our initial problem. We omit here the details and we will next focus on the uniqueness part. 
Set 

Xt := {L^H^ n LlH^ f] LI^bI^^) x L^H-\ for some p G [3, oo[. 

Let (f *, p*) G AV, 1 < 2 < 2 be two solutions of the system ([1]) with the same initial 
data (fo, Oq) and denote 5v = v"^ — f ^, 56 = 6"^ — 6^. Then 

{dt5v - A6v = -P{v^ ■ V5v) - V{5v ■ Wv^) + V{5pe,) 
dt6p + v^ ■ V6p - kA6p= -6v ■ Vp^ 
divf* = 0. 
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Using the maximal smoothing of the heat operator combined with Holder inequality we get 

(38) +\\5p\\lIh-^. 
We have used the classical law product 

\\f9\\L^<\\f\\H49\\^h- 
To estimate we will use Proposition 3.1 of jSj : for every p G [2, oo[ 

< C\\5vVp'\\,.H^.expiC\\Vv'\\ ^ 

We remark that the proof of this result was done in the inviscid case but it can be extended 

to the viscous case with uniform bounds with respect to the parameter k. 

As div 6v = 0, then using Holder inequality and Sobolev embedding ^ yield 

< IIpo||l3||5^^||lihi. 

Thus we get 

< C||po||L3exp(C||Vt;2|| ^ i )\\6v\\^.^.. 
By plugging this estimate into fl38|) . we finally get 

IMl^m ^ {[ \\v\r)\\U\Hr)\\%dr + £ \\6v{T)\\%\\v\r)r^,dr 
+ ||po|li3exp(C||Vi;2|| ^ |)t f \\5v{r)\\l,dT. 

LlBl, Jo 

Since f ^ G LfL°° and G LfH^ then we get the uniqueness by using Gronwall inequality. 
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